Transcription-induced domains form the elementary constraining building blocks of bacterial chromosomes.

Transcription generates local topological and mechanical constraints on the DNA fiber, leading to the generation of supercoiled chromosome domains in bacteria. However, the global impact of transcription on chromosome organization remains elusive, as the scale of genes and operons in bacteria remains well below the resolution of chromosomal contact maps generated using Hi-C (~5-10 kb). Here we combined sub-kb Hi-C contact maps and chromosome engineering to visualize individual transcriptional units. We show that transcriptional units form discrete three-dimensional transcription-induced domains that impose mechanical and topological constraints on their neighboring sequences at larger scales, modifying their localization and dynamics. These results show that transcriptional domains constitute primary building blocks of bacterial chromosome folding and locally impose structural and dynamic constraints.

Design, construction, and functional characterization of a tRNA neochromosome in yeast.

Here, we report the design, construction, and characterization of a tRNA neochromosome, a designer chromosome that functions as an additional, de novo counterpart to the native complement of Saccharomyces cerevisiae. Intending to address one of the central design principles of the Sc2.0 project, the ∼190-kb tRNA neochromosome houses all 275 relocated nuclear tRNA genes. To maximize stability, the design incorporates orthogonal genetic elements from non-S. cerevisiae yeast species. Furthermore, the presence of 283 rox recombination sites enables an orthogonal tRNA SCRaMbLE system. Following construction in yeast, we obtained evidence of a potent selective force, manifesting as a spontaneous doubling in cell ploidy. Furthermore, tRNA sequencing, transcriptomics, proteomics, nucleosome mapping, replication profiling, FISH, and Hi-C were undertaken to investigate questions of tRNA neochromosome behavior and function. Its construction demonstrates the remarkable tractability of the yeast model and opens up opportunities to directly test hypotheses surrounding these essential non-coding RNAs.

Synthetic chromosome fusion: Effects on mitotic and meiotic genome structure and function.

We designed and synthesized synI, which is ∼21.6% shorter than native chrI, the smallest chromosome in Saccharomyces cerevisiae. SynI was designed for attachment to another synthetic chromosome due to concerns surrounding potential instability and karyotype imbalance and is now attached to synIII, yielding the first synthetic yeast fusion chromosome. Additional fusion chromosomes were constructed to study nuclear function. ChrIII-I and chrIX-III-I fusion chromosomes have twisted structures, which depend on silencing protein Sir3. As a smaller chromosome, chrI also faces special challenges in assuring meiotic crossovers required for efficient homolog disjunction. Centromere deletions into fusion chromosomes revealed opposing effects of core centromeres and pericentromeres in modulating deposition of the crossover-promoting protein Red1. These effects extend over 100 kb and promote disproportionate Red1 enrichment, and thus crossover potential, on small chromosomes like chrI. These findings reveal the power of synthetic genomics to uncover new biology and deconvolute complex biological systems.

Karyotype engineering reveals spatio-temporal control of replication firing and gene contacts.

Eukaryotic genomes vary in terms of size, chromosome number, and genetic complexity. Their temporal organization is complex, reflecting coordination between DNA folding and function. Here, we used fused karyotypes of budding yeast to characterize the effects of chromosome length on nuclear architecture. We found that size-matched megachromosomes expand to occupy a larger fraction of the enlarged nucleus. Hi-C maps reveal changes in the three-dimensional structure corresponding to inactivated centromeres and telomeres. De-clustering of inactive centromeres results in their loss of early replication, highlighting a functional correlation between genome organization and replication timing. Repositioning of former telomere-proximal regions on chromosome arms exposed a subset of contacts between flocculin genes. Chromatin reorganization of megachromosomes during cell division remained unperturbed, and it revealed that centromere-rDNA contacts in anaphase, extending over 0.3 Mb on wild-type chromosome, cannot exceed ∼1.7 Mb. Our results highlight the relevance of engineered karyotypes to unveiling relationships between genome organization and function.

Genomics enters the deep learning era.

The tremendous amount of biological sequence data available, combined with the recent methodological breakthrough in deep learning in domains such as computer vision or natural language processing, is leading today to the transformation of bioinformatics through the emergence of deep genomics, the application of deep learning to genomic sequences. We review here the new applications that the use of deep learning enables in the field, focusing on three aspects: the functional annotation of genomes, the sequence determinants of the genome functions and the possibility to write synthetic genomic sequences.

Nucleosome Positioning on Large Tandem DNA Repeats of the '601' Sequence Engineered in Saccharomyces cerevisiae.

The artificial 601 DNA sequence is often used to constrain the position of nucleosomes on a DNA molecule in vitro. Although the ability of the 147 base pair sequence to precisely position a nucleosome in vitro is well documented, application of this property in vivo has been explored only in a few studies and yielded contradictory conclusions. Our goal in the present study was to test the ability of the 601 sequence to dictate nucleosome positioning in Saccharomyces cerevisiae in the context of a long tandem repeat array inserted in a yeast chromosome. We engineered such arrays with three different repeat size, namely 167, 197 and 237 base pairs. Although our arrays are able to position nucleosomes in vitro, analysis of nucleosome occupancy in vivo revealed that nucleosomes are not preferentially positioned as expected on the 601-core sequence along the repeats and that the measured nucleosome repeat length does not correspond to the one expected by design. Altogether our results demonstrate that the rules defining nucleosome positions on this DNA sequence in vitro are not valid in vivo, at least in this chromosomal context, questioning the relevance of using the 601 sequence in vivo to achieve precise nucleosome positioning on designer synthetic DNA sequences.

Contribution of 3D genome topological domains to genetic risk of cancers: a genome-wide computational study.

BACKGROUND: Genome-wide association studies have identified statistical associations between various diseases, including cancers, and a large number of single-nucleotide polymorphisms (SNPs). However, they provide no direct explanation of the mechanisms underlying the association. Based on the recent discovery that changes in three-dimensional genome organization may have functional consequences on gene regulation favoring diseases, we investigated systematically the genome-wide distribution of disease-associated SNPs with respect to a specific feature of 3D genome organization: topologically associating domains (TADs) and their borders. RESULTS: For each of 449 diseases, we tested whether the associated SNPs are present in TAD borders more often than observed by chance, where chance (i.e., the null model in statistical terms) corresponds to the same number of pointwise loci drawn at random either in the entire genome, or in the entire set of disease-associated SNPs listed in the GWAS catalog. Our analysis shows that a fraction of diseases displays such a preferential localization of their risk loci. Moreover, cancers are relatively more frequent among these diseases, and this predominance is generally enhanced when considering only intergenic SNPs. The structure of SNP-based diseasome networks confirms that localization of risk loci in TAD borders differs between cancers and non-cancer diseases. Furthermore, different TAD border enrichments are observed in embryonic stem cells and differentiated cells, consistent with changes in topological domains along embryogenesis and delineating their contribution to disease risk. CONCLUSIONS: Our results suggest that, for certain diseases, part of the genetic risk lies in a local genetic variation affecting the genome partitioning in topologically insulated domains. Investigating this possible contribution to genetic risk is particularly relevant in cancers. This study thus opens a way of interpreting genome-wide association studies, by distinguishing two types of disease-associated SNPs: one with an effect on an individual gene, the other acting in interplay with 3D genome organization.

The 3D Organization of Chromatin Colors in Mammalian Nuclei.

While many computational methods have been proposed for 3D chromosome reconstruction from chromosomal contact maps, these methods are rarely used for the interpretation of such experimental data, in particular Hi-C data. We posit that this is due to the lack of an easy-to-use implementation of the proposed algorithms, as well as to the important computational cost of most methods. We here give a detailed implementation of the fast ShRec3D algorithm. We provide a tutorial that will enable the reader to reconstruct 3D consensus structures for human chromosomes and to decorate these structures with chromatin epigenetic states. We use this methodology to show that the bivalent chromatin, including Polycomb-rich domains, is spatially segregated and located in between the active and the quiescent chromatin compartments.

Improving distance measures between genomic tracks with mutual proximity.

An increasing number of genomic tracks such as DNA methylation, histone modifications or transcriptomes are being produced to annotate genomes with functional states. The comparison of such high dimensional vectors obtained under various experimental conditions requires the use of a distance or dissimilarity measure. Pearson, Cosine and $L_{p}$-norm distances are commonly used for both count and binary vectors. In this article, we highlight how enhancement methods such as the contrast increasing mutual proximity' (MP) or local scaling' improve common distance measures. We present a systematic approach to evaluate the performance of such enhanced distance measures in terms of separability of groups of experimental replicates to outline their effect. We show that the MP' applied on the various distance measures drastically increases performance. Depending on the type of epigenetic experiment, MP' coupled together with Pearson, Cosine, $L_1$, Yule or Jaccard distances proves to be highly efficient in discriminating epigenomic profiles.

Epigenetic rewriting at centromeric DNA repeats leads to increased chromatin accessibility and chromosomal instability.

BACKGROUND: Centromeric regions of human chromosomes contain large numbers of tandemly repeated α-satellite sequences. These sequences are covered with constitutive heterochromatin which is enriched in trimethylation of histone H3 on lysine 9 (H3K9me3). Although well studied using artificial chromosomes and global perturbations, the contribution of this epigenetic mark to chromatin structure and genome stability remains poorly known in a more natural context. RESULTS: Using transcriptional activator-like effectors (TALEs) fused to a histone lysine demethylase (KDM4B), we were able to reduce the level of H3K9me3 on the α-satellites repeats of human chromosome 7. We show that the removal of H3K9me3 affects chromatin structure by increasing the accessibility of DNA repeats to the TALE protein. Tethering TALE-demethylase to centromeric repeats impairs the recruitment of HP1α and proteins of Chromosomal Passenger Complex (CPC) on this specific centromere without affecting CENP-A loading. Finally, the epigenetic re-writing by the TALE-KDM4B affects specifically the stability of chromosome 7 upon mitosis, highlighting the importance of H3K9me3 in centromere integrity and chromosome stability, mediated by the recruitment of HP1α and the CPC. CONCLUSION: Our cellular model allows to demonstrate the direct role of pericentromeric H3K9me3 epigenetic mark on centromere integrity and function in a natural context and opens interesting possibilities for further studies regarding the role of the H3K9me3 mark.

MNHN-Tree-Tools: a toolbox for tree inference using multi-scale clustering of a set of sequences.

SUMMARY: Genomic sequences are widely used to infer the evolutionary history of a given group of individuals. Many methods have been developed for sequence clustering and tree building. In the early days of genome sequencing, these were often limited to hundreds of sequences but due to the surge of high throughput sequencing, it is now common to have millions of sampled sequences at hand. We introduce MNHN-Tree-Tools, a high performance set of algorithms that builds multi-scale, nested clusters of sequences found in a FASTA file. MNHN-Tree-Tools does not rely on multiple sequence alignment and can thus be used on large datasets to infer a sequence tree. Herein, we outline two applications: a human alpha-satellite repeats classification and a tree of life derivation from 16S/18S rDNA sequences. AVAILABILITY AND IMPLEMENTATION: Open source with a Zlib License via the Git protocol: https://gitlab.in2p3.fr/mnhn-tools/mnhn-tree-tools. MANUAL: A detailed users guide and tutorial: https://gitlab.in2p3.fr/mnhn-tools/mnhn-tree-tools-manual/-/raw/master/manual.pdf. WEBSITE AND FAQ: http://treetools.haschka.net. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

keras_dna: a wrapper for fast implementation of deep learning models in genomics.

SUMMARY: Prediction of genomic annotations from DNA sequences using deep learning is today becoming a flourishing field with many applications. Nevertheless, there are still difficulties in handling data in order to conveniently build and train models dedicated for specific end-user's tasks. keras_dna is designed for an easy implementation of Keras models (TensorFlow high level API) for genomics. It can handle standard bioinformatic files formats as inputs such as bigwig, gff, bed, wig, bedGraph or fasta and returns standardized inputs for model training. keras_dna is designed to implement existing models but also to facilitate the development of news models that can have single or multiple targets or inputs. AVAILABILITY AND IMPLEMENTATION: Freely available with a MIT License using pip install keras_dna or cloning the github repo at https://github.com/etirouthier/keras_dna.git. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genome-wide prediction of DNA mutation effect on nucleosome positions for yeast synthetic genomics.

Genetically modified genomes are often used today in many areas of fundamental and applied research. In many studies, coding or noncoding regions are modified in order to change protein sequences or gene expression levels. Modifying one or several nucleotides in a genome can also lead to unexpected changes in the epigenetic regulation of genes. When designing a synthetic genome with many mutations, it would thus be very informative to be able to predict the effect of these mutations on chromatin. We develop here a deep learning approach that quantifies the effect of every possible single mutation on nucleosome positions on the full Saccharomyces cerevisiae genome. This type of annotation track can be used when designing a modified S. cerevisiae genome. We further highlight how this track can provide new insights on the sequence-dependent mechanisms that drive nucleosomes' positions in vivo.

Genome annotation across species using deep convolutional neural networks.

Application of deep neural network is a rapidly expanding field now reaching many disciplines including genomics. In particular, convolutional neural networks have been exploited for identifying the functional role of short genomic sequences. These approaches rely on gathering large sets of sequences with known functional role, extracting those sequences from whole-genome-annotations. These sets are then split into learning, test and validation sets in order to train the networks. While the obtained networks perform well on validation sets, they often perform poorly when applied on whole genomes in which the ratio of positive over negative examples can be very different than in the training set. We here address this issue by assessing the genome-wide performance of networks trained with sets exhibiting different ratios of positive to negative examples. As a case study, we use sequences encompassing gene starts from the RefGene database as positive examples and random genomic sequences as negative examples. We then demonstrate that models trained using data from one organism can be used to predict gene-start sites in a related species, when using training sets providing good genome-wide performance. This cross-species application of convolutional neural networks provides a new way to annotate any genome from existing high-quality annotations in a related reference species. It also provides a way to determine whether the sequence motifs recognised by chromatin-associated proteins in different species are conserved or not.

The 3D Genome Shapes the Regulatory Code of Developmental Genes.

We revisit the notion of gene regulatory code in embryonic development in the light of recent findings about genome spatial organization. By analogy with the genetic code, we posit that the concept of code can only be used if the corresponding adaptor can clearly be identified. An adaptor is here defined as an intermediary physical entity mediating the correspondence between codewords and objects in a gratuitous and evolvable way. In the context of the gene regulatory code, the encoded objects are the gene expression levels, while the concentrations of specific transcription factors in the cell nucleus provide the codewords. The notion of code is meaningful in the absence of direct physicochemical relationships between the objects and the codewords, when the mediation by an adaptor is required. We propose that a plausible adaptor for this code is the gene domain, that is, the genome segment delimited by topological insulators and comprising the gene and its enhancer regulatory sequences. We review recent evidences, based on genome-wide chromosome conformation capture experiments, showing that preferential contact domains found in metazoan genomes are the physical traces of gene domains. Accordingly, genome 3D folding plays a direct role in shaping the developmental gene regulatory code.

Turing-like patterns in an asymmetric dynamic Ising model.

To investigate novel aspects of pattern formation in spin systems, we use a mapping between reactive concentrations in a reaction-diffusion system and spin orientations in a dynamic multiple-spin Ising model. While pattern formation in Ising models always relies on infinite-range interactions, this mapping allows us to design a finite-range-interactions Ising model that can produce patterns observed in reaction-diffusion systems including Turing patterns with a tunable typical length scale. This model has asymmetric interactions and several spin types coexisting at a site. While we use the example of genetic regulation during embryogenesis to build our model, it can be used to study the behavior of other complex systems of interacting agents.

High-salt-recovered sequences are associated with the active chromosomal compartment and with large ribonucleoprotein complexes including nuclear bodies.

The mammalian cell nucleus contains numerous discrete suborganelles named nuclear bodies. While recruitment of specific genomic regions into these large ribonucleoprotein (RNP) complexes critically contributes to higher-order functional chromatin organization, such regions remain ill-defined. We have developed the high-salt-recovered sequences-sequencing (HRS-seq) method, a straightforward genome-wide approach whereby we isolated and sequenced genomic regions associated with large high-salt insoluble RNP complexes. By using mouse embryonic stem cells (ESCs), we showed that these regions essentially correspond to the most highly expressed genes, and to cis-regulatory sequences like super-enhancers, that belong to the active A chromosomal compartment. They include both cell-type-specific genes, such as pluripotency genes in ESCs, and housekeeping genes associated with nuclear bodies, such as histone and snRNA genes that are central components of Histone Locus Bodies and Cajal bodies. We conclude that HRSs are associated with the active chromosomal compartment and with large RNP complexes including nuclear bodies. Association of such chromosomal regions with nuclear bodies is in agreement with the recently proposed phase separation model for transcription control and might thus play a central role in organizing the active chromosomal compartment in mammals.

Kinetic Signature of Cooperativity in the Irreversible Collapse of a Polymer.

We investigate the kinetics of a polymer collapse due to the formation of irreversible cross-links between its monomers. Using the contact probability P(s) as a scale-dependent order parameter depending on the chemical distance s, our simulations show the emergence of a cooperative pearling instability. Namely, the polymer undergoes a sharp conformational transition to a set of absorbing states characterized by a length scale ξ corresponding to the mean pearl size. This length and the transition time depend on the polymer equilibrium dynamics and the cross-linking rate. We confirm experimentally this transition using a DNA conformation capture experiment in yeast.

Characterizing meiotic chromosomes' structure and pairing using a designer sequence optimized for Hi-C.

In chromosome conformation capture experiments (Hi-C), the accuracy with which contacts are detected varies due to the uneven distribution of restriction sites along genomes. In addition, repeated sequences or homologous regions remain indistinguishable because of the ambiguities they introduce during the alignment of the sequencing reads. We addressed both limitations by designing and engineering 144 kb of a yeast chromosome with regularly spaced restriction sites (Syn-HiC design). In the Syn-HiC region, Hi-C signal-to-noise ratio is enhanced and can be used to measure the shape of an unbiased distribution of contact frequencies, allowing to propose a robust definition of a Hi-C experiment resolution. The redesigned region is also distinguishable from its native homologous counterpart in an otherwise isogenic diploid strain. As a proof of principle, we tracked homologous chromosomes during meiotic prophase in synchronized and pachytene-arrested cells and captured important features of their spatial reorganization, such as chromatin restructuration into arrays of Rec8-delimited loops, centromere declustering, individualization, and pairing. Overall, we illustrate the promises held by redesigning genomic regions to explore complex biological questions.

Multiscale Structuring of the E. coli Chromosome by Nucleoid-Associated and Condensin Proteins.

As in eukaryotes, bacterial genomes are not randomly folded. Bacterial genetic information is generally carried on a circular chromosome with a single origin of replication from which two replication forks proceed bidirectionally toward the opposite terminus region. Here, we investigate the higher-order architecture of the Escherichia coli genome, showing its partition into two structurally distinct entities by a complex and intertwined network of contacts: the replication terminus (ter) region and the rest of the chromosome. Outside of ter, the condensin MukBEF and the ubiquitous nucleoid-associated protein (NAP) HU promote DNA contacts in the megabase range. Within ter, the MatP protein prevents MukBEF activity, and contacts are restricted to ∼280 kb, creating a domain with distinct structural properties. We also show how other NAPs contribute to nucleoid organization, such as H-NS, which restricts short-range interactions. Combined, these results reveal the contributions of major evolutionarily conserved proteins in a bacterial chromosome organization.